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The acceleration theorem for Bloch electrons in a homoge- 
nous external field is usually presented using quasiclassical 
arguments. In quantum mechanical versions the Heisenberg 

equations of motion for an operator k(t) are presented mostly 
without properly defining this operator. This leads to the 
surprising fact that the generally accepted version of the the- 
orem is incorrect for the most natural definition of k. This 
operator is shown not to obey canonical commutation rela- 
tions with the position operator. A similar result is shown for 
the phase operators defined via the Klein factors which take 
care of the change of particle number in the bosonization of 
the field operator in the description of interacting fermions in 
one dimension. The phase operators are also shown not to 
obey canonical commutation relations with the correspond- 
ing particle number operators. Implications of this fact are 
discussed for Tomonaga-Luttinger type models. 



I. INTRODUCTION 

The problem of the electronic motion in a crystal sub- 
ject to a homogeneous electric field is-.treated in virtually 
all textbooks on solid state physicsaa. Using semiclassi- 
cal arguments one arrives at the fc-space form of Newton's 
law for Bloch waves 



n d ± = P 

dt 



(1) 



where F = eE, with E the electric field, is a spatially 
uniform force F. In more advanced texts discussions are 
presented that this law also holds rigorously as a quan- 
tum mechanical lawuu, when k(t) is interpreted as the 



expectation value of an operator k(t) = e 



:Hti 



iHl 



for 



an arbitrary wave packet \ip). In this paper we will point 
out various subtleties in the definition of the operator 

k, which surprisingly are not sufficiently discussed in the 
literature. 

It is well-known that there are certain mathematical 
problems with the fact, that the perturbation — F ■ x is 
an unbounded operatoitm This has led to a long discus- 
sion about the existence oflocalized eigenstates, so-called 
"Wannier-Stark-ladders"LrEl. It has been shown that the 
problem is very different if one restricts the Hilbert space 
to a finite number of Bloch bands or works in the full 
Hilbert spacecl. As the important aspects of the prob- 
lem show up already in one- dimensional systems we re- 
strict our discussion to this case. A little more mathe- 
matical sophistication than usual in solid state physics is 



required to appreciate the importance of a proper defini- 
tion of the operator k. Take for example the single band 
case discussed in section III and let \E ) be a localized 
Wannier-Stark eigenstate of the Hamiltonian. Then ob- 
viously (E \k(t)\E ) = (Eo\k\Eo), i.e. the expectation 
value is time-independent in contradiction to Eq. (1). 
The solution to this problem has some similarity to the 
subtleties of the j^aiionical commutation relations for a 
particle on a ringErEj. But there the wave function itself 
is defined on a compact space, while this is not the case 
for the Bloch electron in a homogenous field where the pe- 
riodicity enters via the /c-space classification of the eigen- 
states in a periodic potential. A similar problem occurs 
in a proper definition of the so-called 'tKJein-factors" U, 
which in the method of "bosonization"E3il3 take proper 
care of the change of the number of particles, when the 
electronic field operator is expressed in terms of boson 
operators. The U are unitary operators and as such can 
be expressed as U — e lM , where M is a self-adjoint op- 
erator. We will show in this paper, that, given U the op- 
erator M remains undetermined to a larger extent than 
one could naively expect. Butrwhatever choice one takes, 
this operator M can be shownllj not to obey the commu- 
tation relation with the particle number which is widely 
assumecO. 

In section II we discuss the problem of an electron in 
a periodic potential and a homogeneous external field. 
First introducing the translation operator T a by one lat- 
tice spacing, different definitions of the operator k are 
given, which all lead to T a = e~ lka , but different forms 
of Newton's law and commutation relations with the po- 
sition operator. The most obvious choice for k does not 
lead to Newton's law in the form of Eq. (I) and allows to 
resolve the contradiction for the case of Wannier-Stark 
states mentioned above. In section III we discuss the 
problem, when the Hilbert space is restricted to a finite 
number of bands. The results of sections II and III can 
directly be used for the discussion of the Klein operators 
in section IV. We finally discuss the implications in the 
rather different situations where the same mathematical 
problem shows up in section V. 



II. A BLOCH ELECTRON IN A HOMOGENOUS 
FIELD 

In this section we discuss the motion of an electron in a 
periodic potential V(x) = V(x+a), where a is the lattice 
constant, subject to an additional homogeneous external 



f 



field F 



this yields the result expected from Eq. (6) 



H 



P 

2m 



+ V(x) - Fx. 



(2) 



In the attempt to derive the /c-space form of Newtonla 
law it is useful to introduce the translation operatortil 
which (actively) shifts the state by one lattice spacing to 
the right {x\T a \^) = (x — a\tp). The commutator of T a 
with the Hamiltonian is given by 



[T a ,H] = [T a , —Fx] — aFT a 



Therefore the Heisenberg equation of motion ihT a (t) = 
[T a ,H](t) has the solution 



f a {t)=f a e 



-iaFt/h 



(4) 



As it is well known from elementary quantum mechan- 
ics books the most obvious way to present the translation 
operator is T a — e ~ la P/ h ; where p is the momentum oper- 
ator. This is not the appropriate choice in the following. 
We nevertheless show which contradictions one seems to 
produce using naive reasoning. With the translation op- 
erator represented in terms of the momentum operator 
Eq. (4) reads 



-iap(t) /h 



-ia(p+Ftl)/h 



(5) 



Obviously it is not allowed to infer p(t) = p + Ftl from 
this equation, as it contradicts for all V ^ the Heisen- 
berg equation of motion 



Pit) 



dV 

dx 



(x(t))+Fl. 



(6) 



We postpone the discussion of what goes wrong with this 
argument and first discuss the case V = 0, i.e. a system 
without the periodic potential ("empty lattice"), where 
no contradiction arises. The spectral decomposition of 
k = p/h reads 



\k)k(k\dk , 



(7) 



e lfc 7v27r, 



where the |fc) are plane wave states (x\k) 
which are normalized as (k\k'} = 8{k — k'). In order to 
find k(t) we do not proceed the simple way using Eq. (6) 
but calculate \k(—t))(k(—t)\, needed later. The operator 
Po.fc = commutes with the kinetic energy and its 

commutator with x follows using (fc|cc| , !/>) = i-M-{k\xjj) as 
For V = this leads to the equation 



[Po,k,x] = ijk p 0,k 
of motion 



ihP ,k(t) 



(8) 



k(t) = 



\k - Ft/h)k(k - Ft/h\dk ='k + Ft/hi 



(9) 



This shows explicitly that for V = just comparing the 
exponents in Eq. (5) provides the correct answer. 

We now address the question of how much an operator 

k can differ from k defined in Eq. (7) and nevertheless 



( 3 ) yield f a 



-ika 



To answer this we consider an arbi- 
trary projection operator P — P 2 . Simply expanding 
the exponential function leads to the identity 

e tap = 1 + (e M - 1) P. (10) 
For a — 2nn, with n integer this simplifies to 



e 27mip = 1. (11) 



If two operators A and B fulfill e lA = e lB multiplication 



with e 27rnlP = i yields e 



e«"e~""" . If P commutes 
with B, one obtains e lA = e < t3+2mt "> , i.e. A and B 
can differ by 27rnP and still fulfill e iA = e l§ . If B is a 
self-adjoint operator with the spectral decomposition 



B 



E 

b 



\b)b(b\ 



(12) 



every restricted sum (integral) P — ^2' b \b)(b\ is a projec- 
tor which commutes with B. This simple observation is 
essential for the problem addressed in this paper. 

As a first application we write k in the integrand in 
Eq. (7) in the form fc = k + —n, with n integer and 

k G [— 7r/a, iv /a]. Then we can rewrite this equation in 
the form 



*= E 

n— — oo 

oo 

E 

1 — — C 

oo 

E 



n— — oo 
oo 



7r / a 

-7r / a 
7r / a 



iz/a 



2 -^n)(k 
a 



2 ^n)k { k- 
a 



a a 



— n\dk 
a 



2tt 



a 



OO 

E 

i— — oo 



■k I a 
— n/a 

2irn 



2n 



n)(k 



2tt 



-n\dk 



p(r. 



(13) 



where now k^ is a bounded operator, and the projection 

i.e. we also have T a = 
e -iak[ ' _ j n p OS ition representation k^ is a nonlocal 
operator 



operators P^> commute with ki a ^ , 



with the solution P$ k (t) 

\k(-t))(k(-t)\ = \k- Ft/h){k 



Pa.k-Ft/h(0), i-e. 
Ft/h\. Using Eq. (7) 



2 



(#IV>> = 



oo 

- y 

2vr ^ 




x | i ke ik( - x - x ">dk\ ip{x')dx' 



— / <5(x — x' — ma) 

1 J-°° m=-oo 

cos [(x — x')ir/a\ asm [(x — x')it/a\ 



E 



(a; — x') 
(-1) 



iam 



tt(x — x') 2 
[^(x — ma) — iJj(x + ma)] 



ip(x')dx' 



(14) 



where we have used Eq. (Al). In the following we sup- 
press the a-dependence of the operator . In order to 
discuss kbit) it is convenient to introduce the operators 



Pk 



E \ k - 



2tt 



n)(k 



2tt 



-n\ = P, 



k+ — m > 



(15) 



where we have indicated that they are periodic in the 
reciprocal lattice (m € Z) and we have 



r/a 

kt = I P k kdk . 

' — 7r/a 



(16) 



If for V = Owe again use \k(-t)) (k(-t)\ = \k-Ft/h){k- 
Ft/Ti\ we obtain Pk(t) — Pk-Ft/h- Because of the peri- 
odicity Pk — P k+ 2 JLm we have Pk(t) = Pk(t + mT) where 
T = 2irh/(aF), and we can restrict the following discus- 
sion to times t € [0, T] 



kbit) 



7r / a 



Pk-Ft/hkdk 



■k /a 
-it/a+Ft/h 



7r j a 



Pk-Ft/h+2n/ a kdk 



7r I a 

-Tr/a+Ft/fi 



Pk-Ft/h.kdk, 



(17) 



where we have used Eq. (15). Now we can substitute 
k 1 = k — F t /h + 2n/a in the first term and k 1 = k — Ft/Ti 
in the second. This yields 



k b {t) 



7r J a 

ik' + Ft/h-2%/a) Pk'dk' 
-Ft/h 

(k' + Ft/h)P k ,dk' (18) 

— 7r/ a 

Apart from the term proportional to 2ix j a we can recom- 



bine the integrals and obtain using f_l a /a P k dk = l 



k b it) = k b + Ft/hi 



2tt 



7r / a 



P k dk. 



(19) 



■K/a-Ft/U 



Comparing with Eq. (9) we see that the expressions of 
kit) and kbit) in terms of k and kb differ by the last term 
on the rhs of Eq. (19). We will discuss its implications 
after we show that in fact Eq. (19) also holds in the pres- 
ence of the periodic potential. We denote the eigenstates 
of H = p 2 /2m + V(x) by \k, a) 



H \k,a) = e k Jk,a) , 



(20) 



where k is the wave vector in the first Brillouin zone and 
a is the band index. When the states are normalized as 
(k, a\k', a') = 5 a . a i6(k—k') they provide a decomposition 
of the unit operator 



E 



7r / a 



77 j a 



a)(fc, a\dk = 1 



(21) 



In the appendix we present the simple proof that the P k 
can also be expressed in terms of the Bloch states 



Pk = 



OO 

E 

n=— oo 



k 



2tt 



n)(k 



2tt 



= J2\k,^){k,a\ (22) 



The Heisenberg equation of motion for P k (t) reads 



ihPk(t) 



Pk, Hq 



(t)-F 



Pk,x 



it) (23) 



If one uses the expression for P k in terms of the Bloch 
states the first term on the rhs is seen to vanish, while the 

expression using the plane waves shows that P k , x 

i-^P k and Eq. (23) simplifies to tiP k (t) + F^P fe (t) = 
This leads to the solution 

Pk(t)=P, 



k-Ft/Ti 



(0), 



(24) 



which expresses the well known fact that a state with a 
single sharp value of k remains a state with a single i 
value of k, independent of any interband transitions 
The further steps for the calculation of kbit) are identi- 
cal to Eqs. (16-19). Therefore Newton's law for Bloch 
electrons in a homogeneous field has the solution kbit) = 
hit + nT) and for t E [0,T] 

2.7T f n / a 

hit) = h + Ft /hi V / \k, a){k, a\dk 

a a Jir/a-Ft/h 

(25) 

This one of the central results of the present paper. The 
last term on the rhs of Eq. (25) is missing in all pre- 
sentations of the "acceleration law" known to me. Its 
appearance can be traced to the fact that the operators 
P k are periodic in the reciprocal lattice. A similar term 
has been discussed previously in connection with a free 
particle on a ringp. 

If one calculates expectation values with a general wave 
packet = J2p I-i^/a a fiik)\ki (3)dk the expectation 



3 



value of (k b (t)) = (k b {t + nT)) for t e [0,T] is given 
by 



(k(t))=(h)+Ft/h-—Y i r /a 



/a-Ft/Ti 



\ap(k)\ 2 dk, 



(26) 



with (kb) = J2f3 j-^/ a \ a P^)\ 2 kdk. For a wave packet 
initially strongly peaked around a wave number fco, the 
expectation value (kbit)) closely follows the saw-tooth 
curve shown in Fig. 1 except at the "times of the Bragg- 
reflections" . This implies that Eq. (1) holds to a good 
approximation except during these reflections, when k 
is identified-. with (kb(t)}. For an arbitrary wave packet 
this is no£3 the case due to the additional term from 
differentiating the last term on the rhs of Eq. (26). This 
can be clearly seen in Fig. 1. 




FIG. 1. Expectation value of the operator kb(t) as a func- 
tion of 2nt/T — taF/h for a normalized Gaussian wave packet 
IV*) — c o S exp(— an 2 )|n), where \n) are the Wannier states 
in a single band model. The values for a are 5 (dotted), 
2 (dashed), 1 (dashed dotted), 0.2 (long dashed) and 0.05 
(solid). The asymptotic saw tooth curve is also shown (dot- 
ted). 

The additional term on the rhs of Eq. (25) allows to 
resolve the contradictions mentioned in the introduction. 
This will be discussed in the next section. It also leads to 
an additional term in the commutator [x, kb] compared 
to [x, k] — il. If we compare the Heisenberg equation of 

motion ikb(t) — [kb, H] (t) — —F kb, x (t) for t = with 

Eq. (25) differentiated with respect to time we readily 
obtain 



x, kb 



111 sr— \ .IT 7T 

I > \-,a){-,a\ 

a £ — ' " n 



(27) 



i.e. x and kb do not obey canonical commutation rela- 
tions. If one goes over from kb to kb by adding 2irinP/a, 
where P is a projection operator commuting— with kb, ad- 
ditional terms appear on the rhs of Eq. (27) .lia A detailed 



discussion of the additional term is given in the next sec- 
tion, where we also present an alternative derivation of 
the result. 



III. DESCRIPTION IN A RESTRICTED 
HILBERT SPACE 

The description of a Bloch electron in the external field 
simplifies considerably when the dynamics of an initial 
state in the form of a wave packet in a single band is 
considered and interband transitions are neglectedEHa. In 
this approximation the electron undergoes "Bloch oscil- 
lations" i.e. (x(t)) is periodic in time. A first step-beyond 
this simple approximation is a two-band modelt3 which 
provides a good description for the dynamics of electrons 
in high quality semiconductor superlattices, with a pair 
of isolated minibandsEHlEiJ. More generally the iV-band 
approximation amounts to the replacement 



N „7r/a 

Hp — > / \k,ai)e k , ai (k, 

i=l J-ir/a 



cti\dk , (28) 



and the restriction of the operator x to the Hilbert space 
spanned by the N bands. In the following it is useful also 
to work with Wannier states 



m, a) 



7r j a 



7r I a 



— ika 



\k, a)dk, 



(29) 



which are localized around the m-th lattice site, i.e. 
(x\m,a) — a a (x — ma). They are orthonormal 
(to, a\n, (3) = SapSmn and therefore can also be used for 
a decomposition of the unit operator in the restricted 
iV-band Hilbert space. Elementary calculation yields for 
the Wannier matrix elements of the position operator 

(to, a\x\n, (3) = maS mn S a p + O^"" 

/CO 
dx a* (x)x (x+(m- n)a) (30) 
-OO 

with (ft^T™)* = r^~ m , which implies for the Bloch state 
matrix elements using (A.l) 

d 

(k, a\x\k , (3) = S a pi—5 27I / a (k - k) + 5 2n/a (k - k')Q a[3 (k), 

(31) 

with n af3 (k) = Eie~ lkal n l a p = [n 0a {k)\* and 5^ /a {-) 
the 27r/a-periodic delta function defined in (Al). In the 
following we suppress the index "27r/a" of the delta func- 
tions as wave numbers are always restricted to the first 
Brillouin zone. 

The translation operator T a when acting on a Wan- 
nier state changes to by one, i.e. T a \m, a) = \m + 1, a). 
Therefore the restriction Ta N ^ on the TV-band space has 



4 



the simple form 

f(N) 



N oo 

E i 

i—l m— — oo 



^ y ] \m + l,ai){m,ai 

2—1 '/ 

From the inverse relation of Eq. (29) 



\k, a) — 



E 



e lkam \m,a) 



(32) 



(33) 



it follows that f^ N) \k,a) = 
Ta can be written as 



-ika 



fc, 



exp (— iak b ) with 



and therefore 



E 

i=l 



7r / a 



7r / a 



\k, cti){k, a.i\kdk 



b 

7T / a 
-7T / a 



Pjfhdk 



(34) 



The only difference to section II is therefore the restricted 
sum over bands. In the equation of motion for \t) one 
again has to determine the commutator with the position 
operator. As now there is no equivalent to Eq. (22) the 
argument presented in appendix A cannot be used. Using 
Eq. (31) it is straightforward to show that 



7— P W 

dk k 



(35) 



by calculating matrix elements in Bloch states of both 
sides of this equation. Then all further arguments of 



section II can be used and one obtains ki (t) 



(*- 



nT) and for t G [0, T] 



(t) 



2vr 



Vl^Ft/ti 

N m/a 



E 



■K/a—Ft/h 



\k, on){k, ai\dk . (36) 



Again the only difference to section II is the restricted 
sum over bands. The main new physics is that for 
the restricted TV-band Hilbertspace localized eigenstates 



|^) = (fW) \ES) 



of if, 



(N) 



Fx<- N \ with n G Z and 



(j, = l,...N are known to existafl The corresponding 
energies E£ form TV "Wannier-Stark" ladders 



EH = E% - naF. 



This property 

-l 



is 



easily 



(37) 

shown using 



the contradiction of (E™\k ( b N) (t)\E%) = (E"\k { b 
with Eq. (1) mentioned in the introduction. The last 
term on the rhs of Eq. (36) allows to avoid the inconsis- 
tency. If one takes the expectation value of Eq. (36) with 
the normalized eigenstate \E%) and differentiates with re- 
spect to time one finds that d{E^\k [ b N) {t)\Efi)/dt = is 
fulfilled if 



Fal. Now we can resolve 



N 

E 

i=l 



a 

2^ 



(38) 



holds. In order to see that this in fact is true one has 
to examine the time-independent Schrodinger equation 
H( N 1\E%) = E%\E£). Using Eq. (31) it reads in Bloch 
representation 

D / — f, \ 

i—{k,a i \K) = -{ " p J (k,CH\Ejt) 



N 



-^n« .(k, aj \E%) 
j=i 



(39) 



With k — * t this looks like a "time" dependent 
Schrodinger equation with a "time" dependent Hermi- 
tian N x N Hamiltonian matrix. The "conservation of 
probability" immediately shows that the lhs of Eq. (38) 
is fc-independent. Integration over the first Brillouin 
zone allows to use the completeness relation. Because 
of (E£\E£) = 1 this immediately proves Eq. (38). In the 
simplest case N = 1 Eq. (39) is a single first order linear 
differential equation and can be trivially integrated and 
allows an janalytical discussion of the eigenstates \E n ) in 
real spacefill 

It is instructive to write the operator kb in Eq. (34) 
also in the Wannier representation. If one performs the 
simple integration as in Eq. (14) in (m,a\k{ N ^\n,0) = 



bapim, a \k b N ^ \n, a) one obtains 



K \n,a) 



As the \(m, a\k{ \n, a)\ decay as l/|m — n\ for |m — 
n | — > oo , the states kb\n, a) lie in the Hilbcrt space 7ijy 
spanned by the Wannier states of the N bands. If one 
would take Eq. (40) as the definition of the operator kb 
it is straightforward to show using Eqs. (33) and (40) 
that 



£7 



(~l) r 



i(m — n)a 



m, a). 



(40) 



k,a) = 2 



(-1) 



1+1 



1=1 



la 



■ sin kla 



f st (k)\k,a), 



(41) 



where f s t(k) is the periodic saw tooth curve which is 
given by f s t{k) = k in the first Brillouin zone. The use 
of Eq. (40) also allows an alternative derivation of the re- 
stricted Hilbert space commutation relation correspond- 
ing to Eq. (27). For that purpose we decompose the 
position operator using Eq. (30) as x^ N ' = x 
with 



(JV) 

d 



N 

EE i m ' a * 



am{m, oti\ 



(42) 



This part of the position operator, which is diagonal in 
the Wannier representation, leads to the nonvanishing 



commutator with k 



(N) 



As 



Xj\m, ouj) — am\m,ai l 



the 



5 



states 5Bj |m, ccj) as well as the states fc^5^|m, a,) 
are normalizable and therefore elements oiTL^. If on the 
other hand fej is applied first, one obtains 



is introduced, because instead of using fermionic occu- 
pation numbers it is possible to work with TV + Mo + 1- 
particle basis states 



x { pk[ N) \n,ai) 



{-iy 



y 



rn 



\m, on) 



(43) 



iw,jv)=n 



(blY 



-\{0},N), 



(46) 



;>o 



These states are not normalizable and therefore (like the 
Bloch states) strictly speaking not elements of the Hilbert 
space Hn- If one nevertheless formally calculates the 



where the operators b[ with / > 1 are defined as 



commutator of £j 



and k[ it is given by 



1 



E 

= -Mn 



(47) 



UN) 



n, on 



E 

|rc.,cO 



(-lj^lm.ai) 



(-ir^(-i) 



(44) In the limit Mq — > oo these operators obey Bose commu- 



\m,ai) 



The first form shows e.g. that the diagonal elements in 
the Wannier representation of the commutator vanish, 
i.e. it cannot be proportional to the unit operator. If Eq. 
(44) is multiplied by (n, aj| from the right and summa- 
tions over n and i are carried out, the expression for the 
commutator reads 



(N) t(N) 
d > K h 



= i ( 1 N - — V^IA; = -,£*,-)(& = -,ai\ 
\ a tt a , 



(45) 

which is the analog of Eq. (27) for the case of the re- 
stricted Hilbert space. As in section II it could have been 
derived also by using the derivative of Eq. (36) at t = 



and using Eq. (31) to see that commutes with k K b 
In the next section we discuss how the above commu- 
tation relation occurs in a completely different context. 
There the analog of the Wannier states is the starting 
point of the description. 



tation relations [bj , bj, 



and 



[h,b\,] = s lv i, 



(48) 



while for finite A/opGqs . (46) and (48) only hold for low 
lying excited statescJ'Ej. If the fermions have spin, all op- 
erators have additional spin labels. For the point about 
"Klein factors" we want to make, it is sufficient to work 
with a single branch of fermions as in Eqs. (46)-(48). 
The Klein factors appear in the attempt to express the 
annihilation operatop.|ei, terms of the Bose operators 
defined in Eq. (47)li3liilH It turns out to he easier to 
work with the auxiliary field operator ?/;(w)lla which in 
the limit Mq — s- oo is defined as 



m= E 



e ilv ci 



(49) 



Its commutation relations with the Bose operators are 

1 



b m ,ip(v) 



1 



e im >(u). 



(50) 



IV. BOSONIZATION AND KLEIN OPERATORS 

In this section we first present a very short introduc- 
tion to the ideas behind bosonization, which is a very 
successful techniouje for treating interacting fermions in 
one dimensior£3TL3. Consider a system of fermions with 
one-particle basis states |^) and the corresponding anni- 
hilation operators c/, where I £ Z is a quantum number 
running from — Mq tAoo. For spinless particles in a Id 
box with hard walls&3 Mq — -j=4. In the treatment for 
interacting fermions LuttingerE3 introduced the purely 
technical device to add "unphysical" one particle states 
with negative I and considered the limit Mq — > oo. In 
order to avoid mathematical subtleties it is more trans- 
parent to introduce the concept of bosonizatipn for finite 
Mq and to perform the limit Mq — > oo lateni3. The state 
where the lowest one particle states from —Mq to TV are 
occupied is denoted by |{0},TV), where the symbol {0} 



The essential observation for the "bosonization" of i/j is 
the fact that exponentials of Bose operators obey analo- 
gous commutation relations 



-61 



1 



-bt 



(51) 



6 t e Er =1 ^„ 



If one therefore defines the "Klein operator" 0(v) as fol- 
lows 



ip{v) = 0{v)e 



7» 



(52) 



Eqs. (50) and (51) imply that 0{v) commutes with the 
b n and b\ . As the exponentials in Eq. (52) conserve the 
particle number, the Klein operator O(v) has to lower the 
particle number by one as the field operator on the lhs of 
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the equation. In the following an explicit construction of 
0(v) is presented, which considerably simplifies earlier 
approaches^. The first step is to show that the state 
O(v)|{0}, N) has no overlap to excited states. Using 

({mi},N-l\6(v)\{0},N) = (53) 

n^<{o},iv-i|(feO m! o(«)|{o},iv) 

and the fact that 0(v) commutes with the 6/, the rhs 
of Eq. (53) is seen to vanish unless all to/ are zero, as 
b t \{0},N) = 0. This implies 



As Af corresponds to in section III, the same steps 



O(v)\{0},N) = c N (v)\{0},N-l) 



(54) 



where cn(v) is a c-number. In order to determine cpf(v) 
we calculate ({0},7V- l\i>(v)\{0}N) using Eq. (49) and 
compare with the result when Eqs. (52) and (54) are 
used. This yields cn(v) = e lNv . If Eq. (54) is multi- 
plied by Ui( b h m /Vm?- from the left, the fact that 6{v) 
commutes with the b\ completely determines 0{v) 

6(v)\{ mi },N) = e lNv \{ mi },N - 1) . 



(55) 



Using the completeness of the states \{nii},N) one can 
write Eq. (55) in the form 



O(v) = Ue 



\Sfv 



(56) 



where AT is the particle number operator relative to the 
Dirac see 

^=EElW' JV ) jV <W,JV| (57) 

{mi} N 

and the particle number changing part U which is inde- 
pendent of v is given by 



U=J2 \{mi},N-l){{mi},N\ 
{™i} 



(58) 



If one identifies the particle number N with the site index 
and the {to;} with the at in Eq. (32), the operator U 
has the form of a translation operator T_ a , which shifts 
the system to the left by one lattice site. Therefore the 
eigenstates of U are the corresponding "Bloch states" 



\{mi},k) 



E < 

JV=— oo 



ikN 



\{ mi },N). (59) 



Now one can again write U — e lk with 

k = / \{mi}, k)k({mi}, k\dk. 

r t J — 7T 
{mi} 



(60) 



as in calculation of 



lead tc 



[AT,k] = i ( 1-2tt ^2 \{mi},k = Tt)({mi},k = Tr\ , 

{mi} J 

(61) 



which differs from what, is usually assumedOEj by the 
second term on the rhscil. 



V. IMPLICATIONS 

Here we want to discuss if the subtleties in the definitions 
of the k operators and their non-canonical commutation 
relations have practical implications. We begin with the 
case of a Bloch electron in a homogenous field. Let us 
first summarize the main results presented in sections II 
and III. The well known fact that an initial state in the 
form of a linear combination of Bloch states with different 
band indices but the same value of k remains in such a 
form with a sharp value of k, was expressed as Pk(t) = 
Pk-Ft/h{fy (Eq. (24)). But in contrast to common belief 
this does not lead to Newton's law in the form of Eq. (1) 
if k(t) is interpreted as the expectation value of the linear 
operator kbit) with fcf, defined in Eq. (34). The property 
usually not correctly taken care of is the periodicity in 
the reciprocal lattice of the operator Pf. (Eqs. (22) or 
(34)) . An obvious way to avoid contradictions is not 
to work with kb at all, but take Eq. (24) as the basic 
result and use the translation operator T a only, e.g. in 
the form of Eq. (32) without introducing a k operator. 
We briefly show that this is a useful way to proceed in the 
single band case. The Hamiltonian Hq can be expressed 
in terms of fa 



Hn = 



(62) 



where the e n fi are the hopping matrix elements. In the 
following we drop the superscript (1) at the single band 
translation operator f^ ■ Its representation in Eq. (32) 
easily leads to the commutation relation [x, T™] = naT" , 

where x = x} . Therefore the Heisenberg equation of 
motion for £(t) reads 



n>0 



(63) 



With Eq. (4) its solution can be found by integrating 
from to t 



x(t) -x = J2 [e«,oT? (e- 

n>0 



■inaFt/h 



1) /F + h.c. 



(64) 



7 



If one takes an arbitrary wave packet as an initial state 
the expectation value of 5x(t) = x(t) — x shows periodic 
behaviour with time T = 2-Kh/aF, the well known Block 
oscillations. For the case of nearest neighbour hopping 
the expectation value of T a determines the amplitude of 
the oscillation. If the wave packet is well localized in k- 
space the amplitude is given by B/2F, where B is the 
band width, while a Wannier-Stark eigenstate \E^) or a 
Wannier state |m) as a initial state lead to zero amplitude 
of the Bloch oscillation. For the former all expectations 
are time independent, while the Wannier state behaves 
similar to a squeezed state with a periodic modulation of 
((5x(t)) 2 ). Note that for this description of the N = 1 
case no /c-operator was needed. For N > 1 the Bloch 
oscillations generically are only quasiperiodic and nmapt- 
ical methods have to be used to determine their formEj'Eil. 
Again one can completely avoid to introduce k\>. 
In the case of the Klein operators there exist examples 
where the analogous id^a to completely avoid phase op- 
erators was advocatedEJ. But the usual procedure is to 
introduce the phase operators and to neglect the last term 
on the rhs of Eq. (61), i.e. to assume canonical com- 
mutation-relations for the particle number and the phase 
operatorEj in qtder to be able to apply powerful field the- 
ory techniques^. It is known that this procedure leads 
to the correct anticommutation relations for the field op- 
erators Eq. (49),L3 but it should be obvious that it can 
also lead to errors. As an example we discuss dynamical 
properties of a system described by a Hamiltonian 

H = H B + /(AO, (65) 

where Hb is a function of the boson operators b n and 
fejj and / is an arbitrary function specified later. Note 
that the two terms in the Hamiltonian commute. A typ- 
ical example for this form is the Tomonaga-Luttinger 
modeO'olij. As a first exercise we calculate U{t) = 
e iHtjj e -iHt (p U tting h = 1) in two different ways. This 
operator is an important ingredient for the calculation of 
fermionic correlation functions of the model. Using the 
definition Eq. (58) the exact evaluation is simple as the 
|{mz}, N) are eigenstates of N, and Hb commutes with 
U. This yields 

U{t) = fre-W^-ZCtf-U]*. (66) 

Next we calculate U(t) = e lk ^ by calculating k (t) us- 
ing the Heisenberg equation of motion and the wrong 
commutation relation [Af, k] — il, which gives [k , H] = 
—if'(J\f). As N is a constant of motion this yields 
k(t)cR = k — if'(J\f)t, where the index CR is a reminder 
that we have used the wrong commutation relation. This 
shows that the expectation value of k(t)cR increases lin- 
early with time in eigenstates \E® , N) of H, while the 
expectation value of the exact k(t) is time independent. 
If we put the incorrect result k{t)cR into e lk ^ we obtain 



U(t) C R = e^-f'Wt] 

= e lk f a (e- u ti /'(#(«*))*^ 5 ( 67 ) 

where J\f(a) = e iak J\fe- iak and T a is the "time" (a) or- 
dering symbol. Again using the wrong commutation re- 
lation yields dJ\f(a)/da = —il, i.e. N (a) = N — al. As 
the Af(a) for different arguments commute, the ordering 
symbol T a in Eq. (67) can be deleted and performing 
the integration shows that U(t)cR agrees with the exact 
result for U(t) in Eq. (66). For the Tjmonaga-Luttinger 
model / is a second order polynomia£3iiJ'El Then f'(Af) 
is linear in Jv and the Baker-Hausdorff relation can be 
used for an alternative proof of U (t) = U(t)cR. for Hamil- 
tonians of the type presented in Eq. (65). The fact that 
the incorrect procedure gives the exact result for U(t), 
shows that for the calculation of fermionic correlation 
functions one is lucky for this "fixed point model" of in- 
teracting fermions in one dimension. If backscattering 
or impurity scattering are included the Hamiltonian has 
additional terms which are not of the simple form of Eq. 
(65). If one wants e.g. to evaluate the partition function 
for a Hamiltonian H+XH' the expression in perturbation 
theory to all orders in A involves correlation functions for 
the system described by H defined in Eq. (65). The ar- 
guments just presented show that in this case the use of 
the wrong commutation relation does not lead to incor- 
rect results. 

Nevertheless the example presented shows that the 
widely accepted use of the wrong commutation relation 
between the phase and the particle number operator re- 
quires proper care. 
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APPENDIX A: 

In this appendix we present the simple proof of Eq. (22) 
which states that the operators originally defined in 
terms of plane wave states also have a simple represen- 
tation in terms of Bloch states of an arbitrary periodic 
potential. In the proof we use the well-known identity 

oo oo 

e mu = 2n S{u + 2ttI) = 2ttS 27I {u), (Al) 

m— — oo I— — oo 

where the index "27r" indicates the 27r-periodicity of the 
argument of the delta function. We use this relation to 
show how Pk acts in the position representation 
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/OO 1 JXJ 

oo 

= rj- 11, eikam ">p( x - am )- ( A2 ) 

m— — oo 

For a Bloch state \k' , a) Bloch's theorem (x— am\k', a) — 
e- ik ' am {x\k',a) inEq. (A2) yields, again using Eq. (Al) 

(x\P k \k',a) =S 2v/a (k-k')(x\k',a). (A3) 

If we define the operator P k = \k, P)(k, (3\ the matrix 
elements (x\P k \k' ,a) follow immediately as 

{x\P k \k',a)=8 2n/a {k-k'){x\k',a). (A4) 

As (x\ and \k',a) in Eqs. (A3) and (A4) are arbitrary 
this proves P k — Pk, he. Eq. (22). 



As a simple example we take k b = kb + P k dk. 

Using [x,P k ] — —idPk/dk this yields [x,k b ] — 

i [l — (2n/a)(P^/ a + Pk 2 — -Pfci)] • F° r a finite number of 

additional integrals in the definition of fcf, one obtains a 
corresponding additional sum in the commutator. 
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